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ly-j . When the Channel State Information (CSI) is known by both the transmitter and the receiver, beamforming 
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techniques employing Singular Value Decomposition (SVD) are commonly used in Multiple-Input Multiple-Output 
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channel coding is added, both of them can be achieved as long as the code rate R c and the number of employed 

£> ■ 

. £h , subchannels S satisfy the condition R C S < 1. By adding a properly designed constellation precoder, both full 

X- 

diversity and full multiplexing can be achieved for both uncoded and coded systems with the trade-off of a higher 
decoding complexity, e.g., Fully Precoded Multiple Beamforming (FPMB) and Bit-Interleaved Coded Multiple 
Beamforming with Full Precoding (BICMB-FP) without the condition R C S < 1. Recently discovered Perfect 
Space-Time Block Code (PSTBC) is a full-rate full-diversity space-time code, which achieves efficient shaping and 
high coding gain for MIMO systems. In this paper, a new technique, Bit-Interleaved Coded Multiple Beamforming 
with Perfect Coding (BICMB-PC), is introduced. BICMB-PC transmits PSTBCs through convolutional coded SVD 
systems. Similarly to BICMB-FP, BICMB-PC achieves both full diversity and full multiplexing, and its performance 
is almost the same as BICMB-FP. The advantage of BICMB-PC is that it can provide a much lower decoding 
complexity than BICMB-FP, since the real and imaginary parts of the received signal can be separated for BICMB- 
PC of dimensions 2 and 4, and only the part corresponding to the coded bit is required to acquire one bit metric 
for the Viterbi decoder. 



I. Introduction 

When Channel State Information (CSI) is available at both the transmitter and the receiver, beamforming 
techniques exploiting Singular Value Decomposition (SVD) are applied in a Multiple-Input Multiple- 
Output (MIMO) system to achieve spatial multiplexing^ and thereby increase the data rate, or to enhance 
the performance [3]. However, spatial multiplexing without channel coding results in the loss of the full 
diversity order 01. To overcome the diversity degradation, Bit-Interleaved Coded Multiple Beamforming 
(BICMB), which interleaves the codewords through the multiple subchannels with different diversity 
orders, was proposed Q, (61. BICMB can achieve both full diversity and full multiplexing as long as 
the code rate R c and the number of employed subchannels S satisfy the condition R C S < 1 0, 0|. 
In BUl, IfTOl , X-Codes and Y-Codes were introduced to increase the diversity of multiple beamforming 
which transmits multiple streams. These techniques do not guarantee full diversity when the number 
of transmit or receive antennas is larger than 2, and require relatively high precoding complexity. In 
IfTTTl . Ifl2l . [fT3l , Ifl4ll . it was shown that by employing the constellation precoding technique, which has 
very low precoding complexity, full diversity and full multiplexing can be achieved simultaneously for 
both uncoded and convolutional coded SVD systems with the trade-off of a higher decoding complexity. 
Specifically, in the uncoded case, full diversity requires that all streams are precoded, i.e., Fully Precoded 
Multiple Beamforming (FPMB). A similar result was reported in lfi~5l with a technique employing the 
rotated Quadrature Amplitude Modulation (QAM) constellation. On the other hand, for the convolutional 
coded SVD systems without the condition R C S < 1, other than full precoding, i.e., Bit-Interleaved Coded 
Multiple Beamforming with Full Precoding (BICMB-FP), partial precoding, i.e., Bit-Interleaved Coded 
Multiple Beamforming with Partial Precoding (BICMB-PP) could also achieve both full diversity and full 
multiplexing with the properly designed combination of the convolutional code, the bit interleaver, and 
the constellation precoder. 

In |[T6ll . the Perfect Space-Time Block Code (PSTBC) was introduced for dimensions 2, 3, 4, and 
6. PSTBCs have the full rate, full diversity, nonvanishing minimum determinant for increasing spectral 
efficiency, uniform average transmitted energy per antenna, good shaping of the constellation, and high 
coding gain. In ifTTII . PSTBCs were generalized to any dimension. However, it was proved in lfi~8l that 
particular PSTBCs, yielding increased coding gain, only exist in dimensions 2, 3, 4, and 6. Due to the 
advantages of PSTBCs, the Golden Code (GC), which is the best known PSTBC for MIMO systems 

'in this paper, the term "spatial multiplexing" is used to describe the number of spatial subchannels, as in (T). Note that the term is 
different from "spatial multiplexing gain" defined in (2)- 
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with two transmit and two receive antennas ltT9ll , EDI , has been incorporated into the 802. 16e Worldwide 
Interoperability for Microwave Access (WiMAX) standard fl2T|. 

In our previous work ll22l . Perfect Coded Multiple Beamforming (PCMB) was proposed. PCMB 
combines PSTBCs with uncoded multiple beamforming, and achieves full diversity, full multiplexing, 
and full rate at the same time, in a similar fashion to a MIMO system employing PSTBC and FPMB. 
It was shown that for dimensions 2 and 4, all these three techniques have close Bit Error Rate (BER) 
performance, while the worst-case decoding complexity of PCMB is significantly less than a MIMO 
system employing PSTBC for both low and high Signal-to-Noise Ratio (SNR), and is much lower than 
FPMB for low SNR, which provides the advantage of PCMB. 

In this paper, a new technique with both full diversity and full multiplexing for convolutional coded SVD 
systems, Bit-Interleaved Coded Multiple Beamforming with Perfect Coding (BICMB-PC), is proposed. 
BICMB-PC transmits bit-interleaved codewords of PSTBC, instead of PSTBC codewords without channel 
coding for PCMB, through the multiple subchannels. Diversity analysis of BICMB-PC is carried out to 
prove that it achieves the full diversity order. Simulation results show that BICMB-PC achieves almost the 
same BER performance as BICMB-FP, which is also a technique for convolutional coded SVD systems 
with both full diversity and full multiplexing. Moreover, the decoding complexity analysis shows the 
advantage of BICMB-PC, which has much lower complexity than BICMB-FP for both low and high 
SNR in dimensions 2 and 4. The reason is that the real and imaginary parts of the received signal of 
BICMB-PC can be separated, which is not applied for BICMB-FP, and only the part corresponding to the 
coded bit is required to calculate one bit metric for the Viterbi decoder. Compared to the uncoded system, 
the complexity reduction from BICMB-FP to BICMB-PC is greater than the reduction from FPMB to 
PCMB achieved in E21 for dimensions 2 and 4. Moreover, since the precoded part of BICMB-PP could 
be considered as a smaller dimensional BICMB-FP, BICMB-PC of dimensions 2 and 4 could be applied 
to replace the precoded part and reduce the complexity for BICMB-PP. 

The remainder of this paper is organized as follows: In Section HH the description of BICMB-PC 
is given. In Section Unl the diversity analysis of BICMB-PC is provided. In Section [IVj the decoding 
technique and complexity analysis of BICMB-PC are shown. In Section |Vj simulation results are shown. 
Finally, a conclusion is provided in Section |VIJ 
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II. BICMB-PC Overview 

The structure of BICMB-PC is presented in Fig. \T\ First, the convolutional encoder of code rate R c , 
possibly combined with a perforation matrix ll23l for a high rate punctured code, generates the bit codeword 
c from the information bits. Then, a random bit-interleaver is applied to generate the interleaved bit 
sequence, which is then modulated by M-QAM or M -HEX [24] and mapped by Gray encoding. Then D 2 
consecutive complex-valued scalar symbols are encoded into one PSTBC codeword, where D 6 {2, 3, 4, 6} 
is the system dimension. Hence, the A;th PSTBC codeword Z k is constructed as 

D 

Z fc = ^diag(Gx t) , fc )E i '- 1 , (1) 

v=l 

where G is an D x D unitary matrix, x„ ^ is an D x 1 vector whose elements are the vth D input 
modulated scalar symbols, and 

1 ••• 

oo i ■■• ; 

o 1 

g ••• 

i, D = 2,4, 
< e¥, £ = 3, 
^ -e 2 ?, D = 6, 

and diag(w = [wi, . . . ,%] T ) denotes a diagonal matrix with diagonal entries w\, . . . , wjj. The selection 
of the G matrix for different dimensions can be found in lfi~6ll . 

The MIMO channel H G C NrXNt is assumed to be quasi-static, Rayleigh, and flat fading, and known by 
both the transmitter and the receiver, where N r and N t denote the number of receive and transmit antennas 
respectively, and C stands for the set of complex numbers. The beamforming matrices are determined 
by the SVD of the MIMO channel, i.e., H = UAV H , where U and V are unitary matrices, and A is 
a diagonal matrix whose sth diagonal element, A s E IR + , is a singular value of H in decreasing order, 
where IR + denotes the set of positive real numbers. When S streams are transmitted at the same time, 
the first S vectors of U and V are chosen to be used as beamforming matrices at the receiver and the 
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transmitter, respectively. In the case of BICMB-PC, N r = N t = S = D. 
The received signal corresponding to the kih. PSTBC codeword is 

Y fc = AZ fc + N fc , (2) 

where Y fc is an D x D complex- valued matrix, and N fc is the D x D complex-valued additive white 
Gaussian noise matrix whose elements have zero mean and variance iVo = D/SNR. The channel matrix 
H is complex Gaussian with zero mean and unit variance. The total transmitted power is scaled as D in 
order to make the received SNR SNR. 

The location of the coded bit c k < within the PSTBC codeword sequence is denoted as k' —» (k, (m,n), j), 
where k, (m, n), and j are the index of the PSTBC codewords, the symbol position in X fe = [x ljfc , . . . , x Djk ], 
and the bit position on the label of the scalar symbol xi m ^n),ki respectively. Let \ denote the signal set of 
the modulation scheme, and let xl denote a subset of \ whose labels have b E {0, 1} in the jth bit position. 
Define the one-to-one mapping from X fe to Z fc as Z k = M{X fc }. By using the location information and 
the input-output relation in ©, the receiver calculates the Maximum Likelihood (ML) bit metrics for 
Cfc/ = b as 

r Mj (Y t) *)= min ||Y fe -AM{X}|| 2 , (3) 

where is defined as 

^(m,n)j _ | X . X(u ^ )=(m n) e x ^ andi (tt , c) ^ (m .„) G x}- 

Finally, the ML decoder, which uses the soft-input Viterbi decoding ||25|| to find a codeword with the 
minimum sum weight, makes decisions according to the rule given by G6ll as 

c = argmin V r( m <"^'(Y fc , c k ,). (4) 

c *— ' 

k' 

III. Diversity Analysis 

Based on the bit metrics in ©, the instantaneous Pairwise Error Probability (PEP) between the 
transmitted codeword c and the decoded codeword c is 

Pr(c^c I H) = Pr ( V min ||Y fe - AM{X}|| 2 > V min ||Y fc - AM{X}|| 2 | H . (5) 

I YCol" 1 '")'' Yen'™'"'' 3 / 



Let dn denote the Hamming distance between c and c. Since the bit metrics corresponding to the same 
coded bits between the pairwise errors are the same, © is rewritten as 



Pr (c — > c 



H) = Pr(V min || Y k - AM{X}|| 2 > V min \\Y k - AM{X}|| 2 | H ] , (6) 

\k',d H X€r H' fc',d H / 



where dH stands for the summation of the dn values corresponding to the different coded bits between 
the bit codewords. 
Define and X^ as 



I 2 - 



X fc = arg min ||Y fe - AM{X] 

(7) 

X fc = arg min ||Y fe - AM{X}|| 2 , 

where c k > is the complement of c k < in binary. It is easily found that X fc is different from X fc since the 
sets that xt m ,ri) belong to are disjoint, as can be seen from the definition of '. In the same manner, 

it is clear that X k is different from X k . With Z k = M{X fc } and Z k = M{X fc }, © is rewritten as 

Pr(c^c(H) = Prl > ||Y, - AZ,||~' > > ||Y fc - AZJ| 2 I . (8) 



H) = Pr ( ||Y* - AZ fc || 2 > ll Y * ~ AZ A 

\k',d H k',d H 



Based on the fact that ||Y fc — AZ fe || 2 > || Y fc — AZ fc || 2 and the relation in ©, equation © is upper-bounded 
by 



Pr(c ->■ c | H) < Pr I £ > ^ ||A(Z fc - Z fc )|| 2 J 

V k',d H J 



(9) 



where £ = J2k> d H Tr[— iZ k — Z^^A^N^ — N^A(Zfc — Z k )\. Since £ is a zero-mean Gaussian random 
variable with variance 2N J2 k > d H II-MZ& ~~ Z fc )|| 2 , © is replaced by the Q function as 



Pr(c -> c | H) < Q ( \l Efc ' rfJ|A 2 ^ ZfcH | . (10) 
By using the upper bound on the Q function Q(x) < \e~ x2 / 2 , the average PEP can be upper bounded as 



Pr (c -> c) = E [Pr (c -> c | H)] 

J2k',d H ll A ( Z fc ~ Z * 



< E 



2 eXP \ AN n 



(11) 



In 11221 . it was shown that 



|AZ fe || 2 = Tr[Z^A H AZ fc ] 

S D 



u=l v=l 



2 

v,k\ j 



(12) 



where g^ denotes the uth row of G. By replacing Z k in (PT21) by Z k — Z k , (fTTT) is then rewritten as 



Pr (c ->■ c) < £ 




(13) 
where 

p«,/c = 5^ igr(x„, fc - x„, fc )i 2 . (14) 

1>=1 

The upper bound in (TTSl can be further bounded by employing a theorem from [|27l which is given below. 

Theorem. Consider the largest S < min(A^, N r ) eigenvalues pi s of the uncorrelated central N r x N t 
Wishart matrix that are sorted in decreasing order, and a weight vector p = [pi, • • • ,p,s] T with non- 
negative real elements. In the high SNR regime, an upper bound for the expression E[exp(— 7 Ef=i Psps)}, 
which is used in the diversity analysis of a number of MIMO systems, is 



E 



exp -7 > p 8 fj, 



<C&^„7)- ( "'-* tl,( '"-' +1> 



where 7 w SM?, £ w « constant, p m i n = min Pi ^o {pi} i=1 , 5 « ?/ze index to the first non-zero element 
in the weight vector. 

Proof: See B271. □ 
Based on the aforementioned theorem, full diversity is achieved if and only if 5 = 1, which is equivalent 
to pi > 0. Note that p\^ > in (fl4l) because all elements in gf are nonzero lfT6ll . and therefore 5 = 1. 
By applying the Theorem to (TT3T) . an upper bound of PEP is 



Pr (c -> c) < C 



-N r N t 



AD 



(15) 



Hence, BICMB-PC achieves the full diversity order. 



IV. Decoding 

It was shown in 11221 that each element of AZ^ in d2]) is related to only one of the x^fc. Consequently, 
the elements of AZ^. can be divided into D groups, where the v th group contains elements related to x„ fc , 
and v = 1, • • • , D. 

Take GC (D = 2) as an example, 



AZ> 



Aigfxi ifc Aigfx 2iA; 

zA 2 g 2 r X2 i fc X 2 g2^1,k 



(16) 



The input-output relation in © is then decomposed into two equations as 





^(1,1), k 




Aigfx ljfc 






fl,k = 






+ 




Y{2,2),k 




A2g 2 r xi, fc 




N(2,2),k 



Y2,k 



Y, 



(l,2),fc 



Y(2,l 



(2,i),k 



Aigfx 2>fc 
iA 2 g 2 r x 2i fc 







+ 







(17) 



where Yr m , n ),k and N/ m ,n),k denote the (m, n)th element of and respectively. Let iii^ = [iVn i) N( 2: 2),k] T 
and n 2 fe = [Nn 2 ),k, A r ( 21 ) jfc ] T , then (TT7T ) can be further rewritten as 



y ljk = AGx ljfc + n 1)fc , 
y 2jfc = *AGx 2 t + n 2jfc , 



(18) 



where 



1 

i 



A similar procedure can be applied to larger dimensions. Then in general, the received signal, which 
is divided into D parts, can be represented as 

y v ,k = ^.AGx^fc + h Vtk , (19) 

where v — 1, . . . , D and <f> v = diag(0 1)j i, . . . , 4> v ,d) is a diagonal unitary matrix whose elements satisfy 

1, 1 <u< D + l-v, 
g, D + 2 - v < u < D. 
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By using the QR decomposition of AG = QR, where R is an upper triangular matrix, and the matrix 
Q is unitary, and moving $„Q to the left hand, (fT9l) is rewritten as 

f v ,k = Q H ®vfv,k = Rx^ + n„,fc, (20) 
where h Vi k = Q H ^ii V) k- Then the ML bit metrics in © can be simplified as 

r M ' j (Y,, Cfc ,) = min \\y m , k - Rxf, (21) 
where is a subset of x D , defined as 

^ J = {x = [xi • • • x D } T : x d=n G Xfe, andx^ n G x}- 

The simplified ML bit metrics d2B are similar to BICMB-FP presented in JT2]|, OS, lfl4l| . which are used 
to calculate \M D points by exhaustive search for one bit metric. Hence, the complexity is proportional to 
M D , denoted by 0(M D ). Sphere Decoding (SD) is an alternative for ML with reduced complexity [|28l . 
which reduces the average complexity and provides the worst-case complexity of 0(M D ). Moreover, if 
an efficient implementation of a slicer ||29|| is applied, the worst-case complexity is then 0(M D — 1). 

Particularly, it was proved in E2l that R is a real-valued matrix for dimensions 2 and 4, which implies 
that the real and imaginary parts of y m ,k in (|2T|) can be separated, and only the part corresponding to 
the coded bit is required for calculating one bit metric of the Viterbi decoder. As a result, the decoding 
complexity of BICMB-PC can be further reduced. Assume that square M-QAM is used, whose real and 
imaginary parts are Gray coded separately as two \/M-PAM signals. Define 3ft[^" fe f ] an d ^[C^ ] as tne 
signal sets of the real and the imaginary axes of Q'j, respectively. Therefore, the ML bit metrics in (|2TI) 
can be further simplified for dimensions 2 and 4 as 

r (m,n) J(Yfc;Cfc/) = ^ ||^ ] _ R S R [ X ]||2 ) (22 ) 

R[x]eRK%»] 

if the bit position of cy is on the real part, or 

r M ' 3 (Y t) c,) = min |p[y m , fc ] - R3[x]|| 2 , (23) 

9[x]e9K?£] 

if the bit position of cy is on the imaginary part, where 3?[y m ,fc] an d ^[ym,fc] denote the real and imaginary 
parts of y m ,k, respectively. For (|22b and (|23T) . the worst-case decoding complexity is only <D(M% — 0.5) 
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when SD with rounding (or quantization) procedure for the last layer is employed, which is much lower 
than BICMB-FP of 0(M D - 1). 

Note that in ||22l|. for the uncoded SVD systems with both full diversity and full multiplexing, the 
decoding problem of PCMB is to solve two D-dimensional real-valued problems, while that of FPMB is 
an D-dimensional complex-valued problem in the case of dimensions 2 and 4. In the convolutional coded 
case as shown in this paper, the metric calculation problem of BICMB-PC is only one D-dimensional real- 
valued problem compared to one D-dimensional complex-valued problem of BICMB-FP for dimensions 
2 and 4. Therefore, the complexity reduction from BICMB-FP to BICMB-PC is greater than the reduction 
from FPMB to PCMB achieved in [f22| for dimensions 2 and 4. Moreover, since the precoded part of 
BICMB-PP could be considered as a smaller dimensional BICMB-FP, BICMB-PC of dimensions 2 and 
4 could be applied to replace the precoded part and reduce the complexity for BICMB-PP 

Note that (1221) and (1231) could be further simplified with the trade-off of some performance degradation. 
One way is to replace the square of the difference with the absolute value. Another way is to multiply 
R 1 with 3?[y m ,fc] and 3[y m ,A;]> which is actually a modified Zero-Forcing (ZF) method. Since we focus 
on ML decoding, further details of the trade-off between performance and complexity are not provided. 

V. Simulation Results 

Fig. |2] and Fig. [3] show BER-SNR performance comparison of BICMB-PC and BICMB-FP using 
different modulation schemes for R c = 2/3, 2 x 2 systems, and R c = 4/5, 4 x 4 systems, respectively. 
The constellation precoders for FPMB are selected as the best ones introduced in [fTTTl . Simulation results 
show that BICMB-PC and BICMB-FP achieve almost the same performance for both dimensions with 
all the considered modulation schemes. Moreover, the worst-case decoding complexity of (9(Mt — 0.5) 
to get one bit metric for BICMB-PC is much lower than 0{M D - 1) for BICMB-FP. 

In order to measure the decoding complexity, the average number of real multiplications which are the 
most expensive operations in terms of machine cycles, for acquiring one bit metric is calculated at different 
SNR for BICMB-PC and BICMB-FP, respectively. In J30j, J3D, an efficient reduced complexity decoding 
technique was introduced for BICMB-FP, which is applied in this paper. For fair comparisons, a similar 
decoding technique is applied to BICMB-PC. Fig. @] shows the complexity comparisons for BICMB-PC 
and BICMB-FP in dimensions of 2 and 4 using 64-QAM. For the dimension of 2, the complexity of 
BICMB-PC is 0.8 and 0.5 orders of magnitude lower than BICMB-FP at low and high SNR respectively. 
In the dimension of 4 case, the improvements reach 2.2 and 1.6 orders of magnitude. 
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Additionally, Fig. [5] shows the complexity comparisons for uncoded full-diversity full-multiplexing SVD 
systems of FPMB and PCMB in dimensions of 2 and 4 using 64-QAM. In [|30l, [32]|, an efficient reduced 
complexity SD technique was introduced, which is applied in this paper. For the dimension of 2, the 
complexity of GCMB is 0.5 order of magnitude lower than FPMB at low SNR and similar to FPMB 
at high SNR. In the dimension of 4 case, the complexity of GCMB is 1.7 orders of magnitude lower 
than FPMB at low SNR and similar to FPMB at high SNR. Fig. g] and Fig. \5\ show that the complexity 
reduction from BICMB-FP to BICMB-PC is greater than the reduction from FPMB to PCMB achieved 
in 11221 for dimensions 2 and 4. 

VI. Conclusion 

In this paper, BICMB-PC which combines PSTBC and multiple beamforming technique is proposed. 
It is a new technique with both full diversity and full multiplexing for convolutional coded SVD-MIMO 
systems. Diversity analysis and decoding complexity analysis are provided. It is shown that BICMB- 
PC achieves a similar BER performance to BICMB-FP, which is also a full-diversity full-multiplexing 
technique for convolutional coded SVD-MIMO systems. Particularly, for dimensions 2 and 4, because 
only one of the real or imaginary part of the received signal is required to calculate one bit metric for 
the Viterbi decoder, the worst-case decoding complexity of 0(M^ — 0.5) for BICMB-PC is much lower 
than 0(M D - 1) for BICMB-FP, which provides the advantage of BICMB-PC. Compared to the uncoded 
system, the complexity reduction from BICMB-FP to BICMB-PC is greater than the reduction from FPMB 
to PCMB achieved in ll22l for dimensions 2 and 4. Furthermore, since the precoded part of BICMB-PP 
could be considered as a smaller dimensional BICMB-FP, BICMB-PC of dimensions 2 and 4 could be 
applied to replace the precoded part and reduce the complexity for BICMB-PP. 
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Fig. 1. Structure of BICMB-PC. 
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Fig. 2. BER vs. SNR for BICMB-PC and BICMB-FP for R c = 2/3, 2 x 2 systems. 




Fig. 3. BER vs. SNR for BICMB-PC and BICMB-FP for R c = 4/5, 4x4 systems. 
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Fig. 4. Average number of real multiplications vs. SNR for BICMB-PC and BICMB-FP using 64-QAM. 
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Fig. 5. Average number of real multiplications vs. SNR for PCMB and FPMB using 64-QAM. 



